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1. Introduction 

The classical path approximation which is used in 
most line broadening theories, as well as many other 
atomic and molecular problems, is not always clearly 
understood because it is based on rather heuristic 
arguments. Previous derivations have relied more on 
physical insight than mathematical analysis hence 
such features as the neglect of "back reaction" and 
the assumption that the density matrix is "diagonal in 
wave packet states" are not given to concise mathe- 
matical statements. 

In this paper, statistical methods are used in de- 
riving a precise mathematical expression for the 
"thermal average." The approximations used will be 
given as mathematical statements which are easily 
understood both physically and mathematically. It is 
hoped that this approach will be of some pedagogical 
value in understanding the foundations of the usual 
classical path methods. 

In the following paper [1] l , henceforth referred to 
as II, this approach will be extended to show how the 
usual binary collision theories are obtained, and a 
systematic method for improving these theories is 
outlined. 

2 Construction of a Mathematical Model 

The line broadening problem for gases (or plasmas) 
is basically a study of the distribution of spectral radia- 
tion emitted or absorbed by the gas. The emission (or 
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absorption) is usually produced by a small number of 
particles, and the spectral lines are shifted and broad- 
ened by the interaction between these radiators and 
the remainder of the particles in the gas. Since the 
number of radiators is usually a very small fraction of 
the total number of particles in the gas, the average 
distance between radiators will be large enough that 
any direct interactions between them are negligible. 
We will therefore regard the gas as being divided into 
a large number of cells in such a way that each cell 
contains one radiator and a large number of perturbing 
particles. We will assume that the particles at the 
border between adjacent cells are effectively at infinity 
with respect to the radiators; that is, the interaction 
between a radiator and a perturber at the edge of a 
cell is negligibly small. The distribution of particles 
at the border between cells will then be uniform (i.e., 
not altered by the presence of the radiators) and we 
may neglect all interactions between cells provided 
that we maintain this uniform distribution near the 
edges of the cells. 

We have now represented the gas by a collection of 
noninteracting cells. Each cell may be described by a 
Hamiltonian of the form 



H = H a + H p +V, 



(1) 



where H a is an unperturbed Hamiltonian for the radia- 
tor, H p is an unperturbed Hamiltonian for the remain- 
ing particles in the cell (henceforth called the perturb- 
ers), and V represents the interaction between the 
radiator and these perturbers. 

The probability per unit time that a cell, initially in 
the state \m >, will make a spontaneous dipole transi- 
tion to another state, |m>, is [2] 



(4co| tJl /3ftc 3 )|<7n|d|n>| 



(2) 
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where \m > and \n > are H eigenstates belonging to the 
eigenvalues E m and E n , a) mn is the frequency difference 
(E m —E n )lfi, and d is the dipole moment of the cell. 
Since we are interested only in line spectra, we may 
exclude the continuum radiation produced by the 
perturbers (bremsstrahlung, etc.) by excluding the 
dipole moment of the perturbers. That is, d will be 
the dipole moment operator for the radiator, thus ig- 
noring all transitions which do not involve a change in 
the dipole moment of the radiator. The power radiated 
when a cell makes a dipole transition from \m > to 
\n > is just fitomnPmn, and the contribution this makes 
to the total power spectrum emitted by the gas is 

n<OmnPmn8(Q) — Q) mn ) or 



1(a)) =S m »S(o) — (Omn)\ < m\d\n > | 2 p w 



(6) 



P((o)=(^a) 4 l3c 3 )8(co-(o mn )\<m\d\n>\ 2 



(3) 



if we use delta functions as natural lines. If one wishes 
to use very narrow Lorentzian natural lines [3], it is 
necessary only to use a complex frequency variable 
in place of co; this is of little practical importance 
although it is a mathematical convenience for some 
problems. In P((*)), negative values of o (or a) mn ) 
represent absorptive transitions. 

The total power spectrum emitted by the gas is 
obtained by adding the contributions, P(o)) , from each 
of the cells. Since all cells have the same spectrum of 
accessible states [H eigenstates), this will turn out to 
be a weighted sum of eq (3) over all possible initial and 
final states. If the total number of cells (total number of 
radiators in the gas) is Nr and if the probability of 
finding a cell (one radiator plus its perturbers) in the 
initial state \m> is p m , then the number of cells initially 
in the state \m> is NRp m . Since P m ?i, given by eq (2), 
is the probability that a cell initially in the state \m> 
will make a dipole transition to the state \n>, any 
given transition \m>^ \n> will occur in NRp m Pmn of 
the cells in the gas. The power radiated by this group of 
cells is h(x) mn N Rp m Pmn and their contribution to the 
total power spectrum is ho) mn N R p m P myi b(o) — cu w «). We 
obtain the total power spectrum for the gas by summing 
this expression over all possible |ra>— > \n> transi- 
tions. The total power radiated per unit frequency 
interval is thus given by 

P T ((0) = (WRC0 4 l3c 3 )Xmn$((D-(Dmn) \ < m\d\n >\ 2 p m . 

(4) 

This expression is equivalent to summing eq (3) 
over all \m >-* \n > transitions, weighting each term 
by NRp m , the number of cells initially in the state 
\m> . 

In most line broadening problems one is interested 
only in a small range of frequencies which pertains 
to a single broadened line. In such a case the factor 
to 4 in eq (4) does not change appreciably, and most 
of the frequency dependence of P(o)) is contained 
in the sum over d(o) — a>mn). We thus define a line 
shape function, /(co), by 



The evaluation of 1(a)) is the central problem of 
line broadening theories. Since eq (6) is not multi- 
plied by Nr, we see that it is necessary to consider 
only a single cell (with appropriate weighting of the 
initial states) when evaluating 1(a)). For this reason, 
most line broadening theories start out with a mathe- 
matical model which contains a single radiator im- 
mersed in a gas of perturbing particles. We have 
included this section merely to justify the use of such 
a model, and in all following sections we will refer 
to this model rather than the gas as a whole. 

3. The Weak Coupling Approximation 

Using the results of the previous section, it can 
easily be shown [4, 5, 6] that the line shape for 
spectral radiation is given by 



) = (2t7)-i| X 
= 7r- 1 Re J* 



e iMt C(t)dt 



e^C(t)dt 



(7) 



where C(t) is a correlation function which has the 
property C*(t)=C(—t). For dipole radiation, C(t) 
is an autocorrelation function of the dipole operator d 
for the radiator. This function may be given by 



C(t) = Tr{d-T\t)dT(t)p} 



(8) 



where T(t) and p are the time development opera- 
tor [7] and density matrix for the system (a radiator 
and its perturbers) and the trace is over all states 
of the system. 

The autocorrelation function C(t) may be given 
various physical interpretations. Some of these inter- 
pretations are discussed in detail in section 7. 

Following Fano [8, 9], we regard the system as 
being composed of two weakly coupling subsystems; 
subsystem "a" is the radiator (e.g., an atom, molecule, 
etc.) while subsystem "p" is the gas of perturbers. 
The Hamiltonian, //, for the system has been given 
in eq (1) by (H a -\-H p + V) where H a and H p are Hamil- 
tonians describing the radiator and perturber sub- 
systems and V represents a coupling interaction 
between these subsystems [10]. 

We assume that the radiator and the gas of per- 
turbers are statistically independent; that is, the state 
of the perturbers does not depend explicitly on the 
state of the radiator, and vice versa. This is the neglect 
of "back reaction" in the density matrix. A mathemat- 
ical statement of this approximation has been given 
by Fano [9] in terms of the density matrix: 



p = p {a) p ( P\ 



(9) 



P(o))=(WrO) 4 I3c*)I(o)) 



(5) 



where p (r/) and p^ depend only on radiator and per- 
turber variables respectively. As a first approximation 
we assume that each subsystem is in equilibrium: 
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p(a) = p («)(i/ a ) 

p (p)= p (p)(H p ). 



(10) 



For a canonical ensemble, p a exp {—H/kT}: this 
requires that exp{ — V/kT} be replaced by unity 
leaving p oc exp {—(H (l + H p )/kT} which has the 
form required by eqs (9) and (10). The interpretation 
of this approximation is a matter of some importance 
and it will be discussed in detail in section 8. 
I The assumption of statistical independence, eq (9), 
permits one to write the correlation function, eq (8), 
in the form 

C(t)=Tr a {d-Tr p {T\t)dT(t)p^}p^} 

= Tr a {d-<r(t)dT(t) >p^} (11) 

where the subscripts on the trace operations denote 
a trace over states of the indicated subsystem [9]. 
y The bracket (. . .) represents an ensemble aver- 
age over the perturber subsystem. This average 
takes the place of the usual "thermal average." 



4. The Perturber Ensemble 

We have assumed that the radiator and the per- 
turber subsystems are statistically independent. 
The density operator, p^, for the perturber subsys- 
tem will therefore be the same as the density operator 
for an isolated gas of N particles in a volume f" (the 
volume of a cell) maintained at a constant temperature 
T. Given the restrictions (/V, V, T) , there are many 
possible configurations of positions, x=(xi, X2 . . ., 
x/v), and velocities, v= (vi, v 2 , . . ., v.v), for the par- 
ticles which are consistent with these macroscopic 
constraints. Speaking quantum mechanically, we say 
that there are many different state functions, <£>(x, t), 
which could describe such an isolated gas. We may 
index each of these possibilities with a superscript 
(i), and we imagine a collection of perturber gases, 
called an ensemble [11]. Each (£> (0 ( x > t) must be 
a solution of the Schrodinger equation, 



ot 
which is normalized to unity: 



/ 



p«>(x, *) | 2 rfx=l. 



(12) 



(13) 



Due to its mathematical nature, the concept of 
an ensemble is not readily given to simple physical 
interpretations. However, in the following section it 
will be shown that, by neglecting the influence of 
*the radiator on the state of the perturbers which 
surround it (i.e., "back reaction" effects in the wave 
functions), a physical interpretation of the ensemble 
is possible. In that case we may regard each of the 
cells discussed in section 2 as being a member of 
the ensemble (provided that there are enough cells); 



that is, each c^ (,) will then represent one of the pos- 
sible states for the perturbers in a cell and, for a given 
state of the radiator, the superscript i may then be 
used to index the cells. It must be emphasized thai 
in this section we are not representing the wave func- 
tion for the perturbers in the cells by </? (/) . At this point 
we are simply constructing an approximate density 
operator for these perturbers which is based on an 
isolated gas. In the following section we will use 
the (p^ as basis functions in constructing wave func- 
tions for the perturbers in the cells, however, the 
treatment of the wave functions in the following 
section is distinct from the discussion of the statistics 
in this section. 

For an isolated gas, energy would be a constant 
of the motion hence each <£> (,) will be an eigenfunc- 
tion of H p [13]. If we define a complete orthonormal 
set of H p eigenfunctions {/a(x)} by 



H l f k {*)=m k {*) 



i 



f?(*)M*)dx = &*.j, 



(14) 



(15) 



then, from eqs (12) and (13), it is obvious that 

<p«)(x, t)=j)(x) exp (-it&plh). (16) 

Using the notation 

A(x) = <x|A) 

^(x,0 = (xM i] ; /} (17) 

Equation (16) may also be written in the form 

|^(0 ; *>=|;>exp (-itE^lh). (18) 

The probability of finding the gas of perturbers in 
a particular state <^ (/) will be denoted by cf'K and the 
q^ will be normalized so. that 



W -) =i. 



(19) 



The density operator, p (/>) is then defined [9] by a 
sum over all possible states <£> (/> , weighting each by 
its probability of occurrence^/ : 



i 

=2|;>q<o<;|. 



(20) 



For the macroscopic constraints N, Y , and T, the sta- 
tistical weights are given by 

</<» = exp {- E^lkT}! ^ exp {- E^/kT} (21 ) 



where Eip are the H p eigenvalues. Since the states 
tpW are H p eigenfunctions, they are also eigenfunc- 
tions of the operator exp(— // p / kT)ITrp{exp(—H v \ kT)} 



f 5 
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having the eigenvalues # (/) as given by eq (21). Using 
this property and a closure relation for the <£> (?) we 
may combine eqs (20) and (21) to obtain a more 
familiar form for p (/>) . 

p0» = exp (-H P lkT)ITr p {exv (-H p /kT)}- (22) 

It is obvious, from this equation, that the perturber 
density matrix is diagonal between eigenstates of 
H p . Since the states <£> (/) (x, are eigenf unctions of 
H p , the density matrix will be diagonal between 
these states. In the following section it will be assumed 
that (p (/) (x, t) may, under certain conditions, be 
viewed as a product of wave packets for each of the 
perturbers. We will then be able to state that the 
density matrix is "diagonal in wave packet states." 

5. Approximate Wave Functions 

In this section, approximations will be made con- 
cerning the form of the wave functions for the system 
(radiator and perturbers). We will derive the wave 
functions used by Baranger [4], retaining the super- 
script (i) notation of the previous section in order 
to clarify the influence of the perturbers in these 
wave functions. 

The wave functions for the system must be solu- 
tions of the Schrodinger equation 



ih - ¥(R, x, t) = (fl fl + H p + V)V(R, x, t) , 
ot 



(23) 



where x is a 37V-vector (xi, x 2 , . . ., x. v ), which de- 
notes the positions of trie TV perturbers, and R denotes 
the internal coordinates of the radiator (e.g., for an 
atom, R denotes the positions of the orbital elec- 
trons). 

We wish to expand the wave functions ^(R, 
x, t) in terms of some complete set of basis func- 
tions for the system. Since the functions /a (x) , in- 
troduced in the previous section, form a complete 
set for the perturber subsystem and the H a eigen- 
functions, which we shall callgj(R), form a complete 
set for the radiator subsystem, we may expand 
^(R, x, t) in the following manner: 

^(R,x,0=2^a(0&(R).A(x), (24) 

where the a jk (t) are to be determined by using eq (23) 
and some boundary condition on ^P. 

Using eq (16) we may rewrite eq (24) in the form 



•*(R, x, t) = £ x (A ">( R ' 0^ (A) (x, 



(25) 



where the coefficients 



X <«(R, t) =2 ojt(t)«(R) exp (iE\^tjh) (26) 

J 

must be determined by requiring that ^ satisfy the 
Schrodinger equation, eq (23), as well as some bound- 
ary condition. 



At some initial time, f=0, we may require that the 
perturber gas be represented by one of the wave 
functions <p®: 

^(R, x, 0) =x C0 (R, GV'>(x, 0). (27) 

That is, at some instant in time, the coupling inter- 
action V will be small and the perturbers will behave 
like an isolated gas at that instant. Equation (27) 
constitutes the boundary condition for ^. Substi- 
tuting eq (25) into the Schrodinger equation, eq (23), 
and using the orthogonality of the (£> (A '\ we obtain a 
set of coupled equations for the coefficients x (A): 



h %7~= [#« + PWfOlx^ + S 



dt 



where 



F».J)(t> x 



(j) 



j*k 



VW>(t) = (ifft>it\V\<fP>it) 

V^ k Xt)={<pV;t\V\<p«\t)- 



(28) 

(29) 
(30) 



<, 



We will now assume that the state of the perturbers 
is not affected by the interaction V; that is, 



{<p&; t\V\<pQ\t) <{<(&; t\V\v^,t)- 



(31) 



With this approximation the states (p^ l) may be regarded 
as approximate V eigenfunctions, 



y\ <p (i) ;t> = y®\ <p (0 ; t> , 



(32) 



as well as H p eigenfunctions. This relation provides 
a mathematical statement of the neglect of "back 
reaction" in the wave functions. Since the state of 
the perturbers is now unaffected by the presence of 
the atom, the coupling term in eq (28) vanishes and, 
using the boundary condition given in eq (27), we 
readily obtain the approximate wave function 



^(R,x,0 = x a) (R>0^ (0 (x,0, 



where x^ satisfies the equation 






«) = 



[ifo+J^Ok 



(i) 



(33) 



(34) 
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Comparing eqs (12), (29), (33), and (34) with eqs (16), 
(17), (21), and (22) of reference [4], we see that the 
approximate wave function ^ (/) is identical with the 
approximate wave function introduced by Baranger. 
The superscript (i) now indicates the dependence of 
the total wave function, \P^, on the various possible 
perturber configurations in a cell (as mentioned in 
the previous section). Although this dependence is 
implied by Baranger, it is not illustrated by such a 
parameter. 

It is obvious that a function ^ (/) (R, x, t) having 
the form given in eq (33) cannot rigorously satisfy the 
Schrodinger, equation, eq (23), because the latter 



is, in general, not separable equation. It would be 
possible to find solutions of the form \ (?) (R, x, £)^> (/) (x, 
t), but the "radiator wave function," )^\ would then 
depend explicitly on the perturber coordinates, x, 
as well as R and t. By neglecting "back reaction," 
according to eq (32), we have obtained an approxi- 
mate wave function of the form x (i) <P (i) where x (>) 
does not depend explicitly on perturber coordinates. 
The superscript (i) which appears on this "radiator 
wave function" indicates that x (/) does depend im- 
plicitly on the state of the perturbers through the 
effective potential V^(t); that is, we get different 
functions F (/) for each of the possible states <fM\ 

A general procedure for approximating the solution 
of the Schrodinger equation (23) by a function ^^ 
having the form stated in eq (33) would be an iterative 
one. The first approximation uses <p {i) and x (/) as given 
by eqs (12) and (34). If we chose to, we could obtain 
the next approximation to (p (,) by substituting y i r{i) = 
^ (i *X into the Schrodinger equation (23), using eq 
(34) to remove some of the x (i) terms. Multiplying by 
X (/) * and integrating over dR then gives a differential 
equation for ip {i) which is similar in form to eq (34). 
That is, in addition to the H p which appears in eq (12) 
we would have a potential energy function which is 
averaged over R. The <p {i) solution of this new equation 
is presumably an improvement over the solution of eq 
(12) since the former contains some back reaction 
effects; however, the improved ^? (,) would no longer be 
an H p eigenfunction. Since it is desirable that (ft® be 
an H p eigenfunction and thus diagonalize p (p) , we will 
use the first approximation to <£> (/) (i.e., the solution of 
eq (12)). If back reaction effects are important, it is 
not sufficient to use this first approximation. If a higher 
order approximation is used, the density matrix, p {p \ 
will not be diagonal, as noted by Griem [14]. [In sec. 8 
it will be shown that when these back reaction effects 
become important, the classical path approximation is 
breaking down anyway.] In this case, it would be more 
convenient to use the H p eigenfunctions as a basis set 
in evaluating the trace over perturber coordinates 
(i.e., the thermal average as defined by eq (11)). The 
perturber density matrix would then be diagonal, but 
it would be necessary to make a fully quantum me- 
chanical calculation of the time development operator 
(similar to that discussed by Baranger [4, 5]). Such a 
procedure would not be a classical path calculation 
and we will not consider this case any further. 

In keeping with the usual classical path treatment, 
we next assume that each wave function <^ (,) (x, t) 
may be expressed as a product of wave packet states 
(p^ixj, t) for each of the perturbers: 



having eigenvalues E$ and V (i) (t) which are classical 
functions of perturber coordinates. The validity of this 
approximation is discussed in section (8). 

We now note that the influence of the perturbers 
which is implicit in V {i) {t) may be specified by giving 
the locations, xj'\ of the peaks of the packets <py\ as 
well as their group velocities, vj°. We thus complete 
the classical path picture by regarding the wave 
packets y^ as classical particles described by their 
positions 

x (/) = (x^, x£>, . . ., x$) and velocities 

v<*>= (vi f \ vf>, . . ., v$). 

Since the wave packets (pW are assumed to be sharply 
localized, we may replace V (i) (t) by a classical poten- 
tial function. That is, 



{ 



*>«*(x, 0*V°(x, t) dx= P>(R, x«>, v<'>, t) 



(36) 



where V (i) (R, x (/) , v (/) , t) is a classical potential func- 
tion of the 3/V-vectors x (/) and v (/) , the operator R, and 
the time t. 

6. The Thermal Average 

To obtain an expression for the "thermal average" 
in eq (11), we must first obtain an expression for the 
time development operator, T(t), for the total system. 
Since we are approximating the wave function 
*(R, x, t) by *<'">(R, x, *), we will need the cor- 
responding time development operator 7^(0 which is 
defined by 



^>(R, x, t) = T*>(t)V®(R 9 x, 0). 



(37) 



From eqs (33) and (34), we see that this operator may 
be written in the form 



TV(t) = exp{-itH p /h} TQ($) 
where T$(t) is the solution of 



(38) 



ih - Jf(t) = [H a + ^(R, x«>, v<», r)]jy (t) (39) 



with the boundary conditions T$(0) = 1. 

Since d (see eq (11)) is the dipole operator for the 
radiator, it will commute with the perturber Hamil- 
tonian H p and we may write 



<p»(x f o=n*w x *'>- 



(35) 



The wave packets, tffp, are assumed to be sharply 
localized during the times of interest in line broaden- 
ing. Although one does not usually think of energy 
eigenstates as being wave packets, it is nonetheless 
possible for a wave packet state to be an approximate 
eigenfunction of both H p and V (as required by eq (32)) 



T® f (t)dTW(t) 



W'{t)*TV(t). (40) 

thermal average" 



Substituting this identity into the 
in eq(ll), we obtain 



(. . .) = Tr p {W i (t)dnP(t) P ^} 

A-A-' 



(41) 
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where \k > denotes an arbitrary basis vector for the 
perturber subsystem. Since T^\t) has no explicit 
dependence on perturber coordinates (i.e., it depends 
on perturbers implicitly through the coordinates 

x (/) , v (?) of the wave packets <£> (/) ), we have 



<*in 01 (0d7l°(0l*'> = 8*.A-7\? t (*)dn (0 



(42) 



Substituting this identity into eq (41), we see that 
we need only the diagonal matrix elements of p (;;) 
in the states \k >. Using the matrix elements of p (/>) 
as given by eq (20) and noting the normalization con- 
dition of eq (16), we obtain 



<. . .> = -2q(')TV(t)dr%t). 



(43) 



To complete our evaluation of the "thermal aver- 
age," we must express this weighted average in a form 
which is amenable to practical calculations. 

Since g (/) is a function of Ety, defined by eq (21), we 
may express ql® as a function of x (,) and v (,) by using 

< <?<*>; t\H p \<pM; t> = Ef(x«\ v U) ) 

j 

where V p denotes the interaction between perturbers. 
If the radiator is an ion in a plasma, we may add a 
Coulomb interaction, F ( >(x (/) ), to H p [10]. This ex- 
pression is consistent with the wave packet interpre- 
tation of <£> 0) (x, 0- 

In the classical limit, the sum over (i) is replaced by 
an integral over x (0 and v (/) . We may thus write 
g ( m the form 

g (0( x (0, v (») = exp {-£(0( x (0 9 vd))lkT}IQ, (45) 
where Q is the classical partition function 

Q = [ [ exp {-E<?(x (/) , v^)/AT}^x<'W''>. (46) 
We may now write the "thermal average" in the form 
<. . .> = <?- ] f |dx<W*exp {-£$(*« v®)lkT} 

rg>t( x ^), v<«>, t)dVi}{*s i \ v» o (47) 

where the time development operator T$(xfi\ v (/) , t) 
is obtained by solving eq (39) with the classical poten- 
tial V {i) (R, x (/) , v (/) , t). This expression will be reduced 
to the more familiar averages over collision variables 
(velocity, impact parameter, etc.) in II. 

7. The Autocorrelation Function 

7.1. The Influence of the Density Matrix and the 
Wave Functions 

Before discussing the validity of the approximations 
concerning the density matrix and the wave functions 



(sec. 8), it is useful to first consider the roles they 
play in determining the line shape. Since the line 
shape has been given as a Fourier transform of the 
autocorrelation function, eq (7), it is convenient to 
study this function as well. 

The autocorrelation function was defined in eq (8) 
by C(t) = Tr{d-d(t)p} where d(t) = T f (t)dT(t) 
and p is an equilibrium density matrix for a radiator 
and its perturbers. Since both p and T(t) are func- 
tions of H (see eq (1)), they will commute at all times 
(i.e.,p is stationary) and we maywrite 



Tr {d • d(t)p} = Tr{d(t ) ■ d(t + t )p} 



(48) 



for any time Jo. This relation is characteristic of a 
stationary random process and it indicates that when 
we speak of an initial time, we are referring to any 
arbitrarily chosen instant, not some time infinitely 
far in the past. In writing C(t)= Tr {d • d(t)p} we 
have arbitrarily chosen t = as our initial time. Auto- 
correlation functions [15] such as C(t) fall to zero 
as t increases; this represents the loss of correlation 
between the dipole moment at the time t = Q and 
its value at some later time. 

We may regard the trace, in the defining relation 
for C{t), as a sum over all possible initial configura- 
tions of perturbers and initial states of the radiator; 
each term in this sum is weighted by the matrix 
element of p which gives the probability of finding 
the corresponding initial state. We may also interpret 
this trace as a sum over cells, with each initial state 
corresponding to a given radiator-perturber con- 
figuration as discussed in section 2. For each initial 
state there will be a dipole moment, d, and at some 
time t later this moment will have changed to d{t). 
Once the initial state has been specified, the temporal 
change in d is determined solely by the time develop- 
ment operator T(t), and is not influenced by the 
statistics (density matrix) in any way. 

We may now think of the line shape as being ob- 
tained by taking the Fourier transform of the d • d(t) 
contribution from each initial state, and then summing 
these Fourier transforms, weighting each term by 
the appropriate density matrix element (the proba- 
bility of occurrence of a given configuration in a cell). 
It is thus apparent that the density matrix can influence 
the line shape only by its weighting of each d • d(t) 
contribution; the wave functions influence I(co) 
only by their effect on the calculation of each d • d(t) 
term. The influence of the density matrix is quite 
distinct from the influence of the wave functions, 
and similarly, approximations relevant to the statistics 
(density matrix) are completely independent of the 
approximations relevant to the dynamics (time develop- 
ment operator or wave functions). 

7.2. Classical Correspondence 

a. Classical Theories 

In discussing any quantum mechanical theory it is 
frequently beneficial to have a classical correspond- 
ence available as an aid to understanding the physical 
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properties underlying the mathematical description. 
In a classical path description of line broadening, the 
perturbers are regarded as classical particles, and 
there is little difficulty in understanding the physical 
properties of this subsystem. The radiator is treated 
quantum mechanically however, and it is useful to 
develop its classical analog to clarify the influence of 
collisions on the spectral radiation. 

In classical line broadening theories one considers 
a classical oscillating dipole which emits a continuous 
wave train (referred to as the light train). The light 
train is described by the classical dipole moment 
d(t), and the line shape is given by a Fourier trans- 
form of the dipole autocorrelation function 
C(t) = <dd(t) > just asineq(7); in this case < . . . > 
denotes a classical average over collisions. In the 
Lorentz interruption theory (see p. 571 of ref. [16]), 
it is assumed that a collision terminates the unper- 
turbed light train d(t) = d exp (— icoot) where coo is 
the unperturbed frequency. The autocorrelation func- 
tion is therefore given by, C(t) = d 2 / P(to)U(t < to) 
exp (— i(x)()t)dto, where P(to) is the probability that a 
collision will occur after a time £o, and U(t< to) is a 
step function which equals 1 before the collision occurs 
(i.e., when t < t {) ) and equals zero after the collision. 
Although the line shape is not always given in this 
form (cf. eq (1.2) of ref. [16] noting that P(to) = v c 
exp (— v c k))), the results are the same. In another 
classical theory Weisskopf, and later Lindholm, as- 
sumed that collisions can change the phase of the 
light train [16, 17]. The line shape is again given by 
the Fourier transform of the light train autocorrelation 
function (see also eq (1.13a) of ref. [16] or eq (11.13) of 
ref. [17]) hence a phase shift decreases the phase 
correlation thereby broadening the Fourier spectrum 
(line shape). In these phase shift theories, d(t ) has the 
form exp [— ia>ot — irj (t ) ] where r)(t) represents the 
phase shifts produced by collisions. The autocorrela- 
tion function is C(t) = d 2 exp (— i(Oot) (exp [—ir)(t)] ) 
where < . . . > denotes an average over phase shifts 
(thermal average). The Weisskopf theory assumed that 
small phase shifts (less than unity) have no effect on 
C(t), whereas large phase shifts (greater than unity) 
destroy the phase correlation of the light train and 
cause C{t) to fall to zero. This effect is only apparent 
after one performs an average over a probability dis- 
tribution of phase shifts. To show this, we consider a 
collision which occurs in the time interval (f, f-f At) 
and produces a phase shift no. If the duration of this 
collision r is very short, r < Af , then we have C(t-\- At) 
= C(t) (exp (— ir)o)) (the two factors are independent 
and are averaged separately as in the impact approxi- 
mation; see sec. (11.3.2) of ref. [17]), where (exp 
( — "?o)> = / exp(—i7]o)P(r} { ))(lr) {) and P( no) denotes the 
probability of finding a phase shift of magnitude no. If 
P(r)o) does not vary rapidly, then the rapid oscillation 
of exp (— irjo) for no > 1 will produce 



f 



exp {—ir)o)P(r)o)d7)o< 1 



(the integral from — w to — 1 is also much less than 1) 



whereas the integral from — 1 to + 1 is on the order of 
unity. Thus, for no > 1, we have C(t-\- At) — and for 
no < 1, C(t + At) —C(t). While the average over no 
is given more accurately by later phase shift theories 
(e.g., Lindholm's theory), the qualitative results for the 
line shape stated by Weisskopf for no > 1 are basically 
unchanged. It is in this sense that we will regard a 
large phase shift as an effective termination of the 
light train. 

b. Quantum Mechanical Theory 

In the quantum mechanical problem the radiators 
in a gas emit a series of spectral lines which are broad- 
ened due to the perturbation V. We are usually interest- 
ed in just one of these broadened lines corresponding 
to some given unperturbed frequency coo. The photons 
which contribute to this particular line are emitted in 
radiative transitions between some particular initial 
and final energy levels of the radiator (which may be 
perturbed). In this case, the light train is the wave 
function for a photon emitted in such a transition; 
the amplitude of the light train gives the probability of 
finding a photon emitted after some time t. Baranger 
has shown that the function C(t) defined in eq (8) may 
be interpreted as the autocorrelation function for the 
light train amplitude (see p. 498 of ref. [4]). Thus the 
line shape is obtained from the Fourier transform of 
the light train autocorrelation function just as in the 
classical theories. 

The effects of a collision may be roughly divided 
into two categories: adiabatic and nonadiabatic. We 
will discuss these separately although a single collision 
may have both effects. An adiabatic collision will 
simply shift the energy levels thereby producing a 
phase shift (which may be time dependent) in the 
time development operator T(t). Such collisions will 
reduce the phase correlation between d(0) and d{t) = 
T f (t)dT(t) just as the collisions in the classical phase 
shift theories. A nonadiabatic collision will induce a 
radiationless transition out of the initial state. This 
type of transition usually "terminates" the probability 
that the radiator will emit a photon with a frequency 
near ct> . However, if there is another state whose energy 
is very close to that of the initial state, nonadiabatic 
collisions may simply produce transitions back and 
forth between these two states without completely 
destroying the probability of emission near coo. The 
latter situation is usually referred to as the problem 
of overlapping lines (although it is sometimes called 
a "phase-memory" effect) [18]. In a classical oscil- 
lator having more than one mode of oscillation, a non- 
adiabatic collision simply changes the oscillation from 
one mode to another. The Lorentz interruption theory 
implicitly assumes that the mode whose frequency is 
coo is the only mode of interest and a (nonadiabatic) 
collision terminates this mode. However, if there were 
another mode with a frequency very close to ojo, non- 
adiabatic collisions could switch the oscillator back 
and forth between these two modes without destroy- 
ing the phase coherence of the light train (overlapping 
line case). For example, an oscillator in the mode o>o 
will have a phase coot at some time t (apart from an 
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arbitrary initial phase). If this oscillator is switched to 
the mode coi at some time t and, after a time £o, is 
switched back to o>o, the resulting phase shift in the 
natural oscillation is r) — [a)o(t-\- to) — (mt-\- (Oito)] = 
(ct)o — (x)i)to. Estimating to by the mean free time {\jv c ) 
and noting that the collision frequency v c may be 
identified with the halfwidth A0J1/2, (in a Lorentz pic- 
ture), we have 77 — (o>o — a>i) /Aa>i/ 2 . If (o>o — o>i ) > Ao>i/2 
(isolated lines), we have 77 > 1 and we may say that 
the light train is terminated by nonadiabatic collisions 
(as the Lorentz theory assumes). If (ojq — wi) < Aa>i/2 
(overlapping lines), we may say that the light train has 
retained a "memory of its initial phase." Thus we 
see that nonadiabatic collisions are analogous to those 
considered by the Lorentz interruption theory only in 
the case of isolated lines. In the case of overlapping 
lines, nonadiabatic collisions are not treated by any 
classical theories although the classical analog is 
obvious. We also note that a collision which produces 
both adiabatic and nonadiabatic effects cannot be 
conveniently treated by classical theories either. 

To make our discussion of the quantum mechanical 
C{t) more quantitative we will consider the case of a 
quantitative three-state radiator having a ground state 
with energy E {) and two excited states with energies E\ 
and Ei. We will also assume that the ground state is not 
perturbed by the interaction V(t) and we use Foo=Foi 
= Fo2 = (as well as Vij=Vjj). A generalization to more 
complicated cases may be achieved by means of an 
obvious extension of the following arguments. The time 
development operator for the radiator (see eq (39)) will 
be written in an interaction representation defined by 



T a (t) = exp(-itHalh)r(t), 



(49) 



where 



T{t) = exp 



- am 



J>(o4 



V (t) = exp (itHalh) V{t) exp (- itHJh) , (50) 

and 6 is the time ordering operator [19]. Notice that 
Vqj=Vjo = implies that 7oo = l and TI)j=Tjq = for 
j # 0; this simplification is the reason for neglecting 
ground state interactions. 

If we are interested only in the line whose unper- 
turbed frequency is (t)\o= (E\—Eo)lh, we do not need 
all of the terms in the trace over states of the radiator 
(see eq (11)). We need only consider the terms 

C(0 = 2 djodoi(T[j)p( a) exp (— ia) U) t) 

j 

+ 2 djod {)2 (T! zj )pW exp (- io* 2 ot) , (51) 



where pj o) a exp (—Ej/kT). The natural oscillation exp 
(— ia>iot) in the first term may be regarded as an un- 
perturbed light train which causes the Fourier trans- 
form of this term (eq (7)) to be large in the vicinity 
of coio. Thus, when calculating the line coio we need 



only consider those terms in eq (11) whose natural 
oscillation is close to o>io. The second term in eq (51) 
describes radiation near CO20; this term will be important 
(when calculating the com line) only if the lines o>io 
and o>20 overlap (i.e., if their halfwidths are on the 
order of (ct>io — CO20)). It will be shown later in this sec- 
tion that the phase of the off-diagonal terms such as 
T' V2 varies according to the line spacing o>i 2 ; thus when 
the average time between collisions (which is on the 
order of I/A001/2) is small compared with l/o>i2, these 
off-diagonal terms will become important and phase 
memory effects (overlapping lines) will have to be 
considered. If coio and CO20 are isolated lines (i.e., their 
halfwidths are much smaller than W12), the second term 
in eq (51) (as well as the off-diagonal T matrix elements) 
may be neglected. For some lines in the microwave 
region it may happen that the halfwidth of o>io is on the 
order of o>io; in this case it would be necessary to add 
the "negative resonance" terms, exp (— icooiO , to 
eq (51). 

If we factor the natural oscillation out of C(0? 
we may write C(t) = exp (— icoiot)C{t) where C(t) 
is a decreasing function of t [15]. The line shape (see 
eq (7)) is proportional to / exp {ithu))C(t)dt where 
Aco= (00 — o>io). Since exp (it&aj) oscillates rapidly 
for t > I/Aw, the intensity at the point Aco is deter- 
mined primarily by the part of C(t) for which 



0^^^ l/Ao>. 



(52) 



For this reason, 1/Aco is frequently referred to as the 
time of interest in evaluating the intensity at the 
point Aw. 

If the duration of a collision is much longer than the 
time of interest, the perturber may be regarded as 
stationary during this time. This quasi-static approxi- 
mation is frequently employed in the description 
of ions in a plasma [4, 6]. Such static perturbations 
simply add a constant to the oscillation frequencies 
coio, CO20, etc. The influence of dynamic perturbations, 
for which the collision time is less than or the order of 
the time of interest, is obtained by calculating (7^) 
and(r4,->. 

We are interested in the classical correspondence 
for the quantum mechanical C(t) primarily as an aid 
in establishing a criterion for the so-called strong col- 
lisions (this will be done in sec. 7.3). When a single 
collision is capable of reducing C{t) to zero, it is re- 
ferred to as a strong collision. To study this effect 
we need only consider the influence of a single col- 
lision on C(t). We will consider some particular 
type of collision (e.g., given values of impact parameter, 
velocity, etc.) which is described by the interaction 
potential V(t). We will assume that this is the first 
collision to occur after some arbitrarily chosen start- 
ing time £ = (if a collision does not reduce C(t) 
to zero, it is "weak," and to find a strong collision 
criterion we may neglect the effects of all previous 
"weak" collisions). The collision will begin at some 
time £0 ^ 0, and we define a collision time r such that 
V(t) = unless to^ t^ to + r. From eq (50) we see 
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that T'(t) = l when t < t {) ; when £>(*o+t), we 
change the variable of integration to t"=(t' — to) 
and obtain 

r(0=exp (itoHJh)S exp (-itoHJh), (53) 

where S is the S-matrix for the collision, 

S=l-(i7ft) I F(*)A 

Jo 

-(l/ft) 2 f T fV(0^')A'*+- ' •, (54) 

Jo Jo 



and 



V(t) = exp (UHalh)V(t + to) exp (-itHJh). (55) 



The contribution to C(0 from this collision will be 
called C(V(t), t); this function will have the same 
form as C(t) in eq (51) except that there will be no 
average < . . . > over collisions. To obtain C(t) 
we would simply average C(V(t), t) over the collision 
parameters for the first collision (and all successive 
collisions if the first collision is not strong). 

c. Isolated Lines 

If the lines a>io and W20 are isolated, C(V(t), t) 
is given by Xjdj do\T' u p$ a) exp (—iaiiot). In this 
series T[ o = (no ground state interaction) and T[ 2 
is unimportant (for isolated lines only). As mentioned 
earlier, to show that T{ 2 is unimportant it is necessary 
to average over to. This has been done for the line 
wings in eq (19) of ref. [20] and the second term in this 
equation, which contains the off-diagonal elements of 
the time development operator, is negligible in the 
case of isolated lines (see also the following section 
on overlapping lines). Using this fact we have 



C(V(t), t)<xdtopiT'n(t)exp(-ia> 10 t). 



(56) 



Since the complex function T n (t) can always be writ- 
ten in the classical form exp [— ir)(t) — e(t)], the 
classical correspondence for C(V(t), t) becomes ob- 
vious. Collisions for which e = are equivalent to the 
collisions in the classical phase shift theories and those 
for which e > 1 are equivalent to the nonadiabatic col- 
lisions in the Lorentz theory. To illustrate the func- 
tional form of 7} and e, we will evaluate T' n (t) , for 
t> (j + t), to lowest order in V(t). We will use 
V\ 1 = ^22 = since V has no diagonal matrix elements 
in most problems (i.e., V usually has non vanishing 
matrix elements only between states with different 
angular momenta). Using eqs (53), (54), and (55) we 
obtain 

r;, (0 = 1 -ft- 2 '£ f Vvz{t+ t0 )v 2 i (t'+to) 

exp {iAE 2 i(*' —t)/h}dt'dt 



» l-i(Ffi/*ABn) P [exp {ttA&u/ft}-l]<&, (57) 

where AE21 = (E 2 — E\) and V 2 \ denotes the average 
value of V 2 \{t) during the collision. For fast collisions 
(A2?2iT/ft) < 1, we have 



T' n (t)~l-HV2lTlh)\ 



(58) 



For slower collisions (A^ir/ft) > 1, we replace r by 
^° when integrating the exponential in eq (57). We 
then replace A£ 2 i by (A^i + id), perform the integral, 
and then take the limit 8— > 0. The justification for this 
procedure is the "radiation damping" which produces 
the natural line width (see sec. 13. A of ref. [3]). In 
this manner we obtain 

r„(0 ~l + i(JXT/ftAfti). (59) 

These results may be summarized by 

C(V(t), t) a dfopi exp (— iaho«+ilJ-€), (60) 

where 

f (AE 2X Tlh)<\ 

7i = \ (61) 

l(^T/ftA£ 2 i) (A£ 21 r/ft) >l 



\HV21Tlh) 2 







(AE 21 T/ft) <^ 1 
(A£«it/*) >1. 



(62) 



Notice that if the perturbation is not strong enough to 
mix adjacent levels, V 2 i < &E 2 \, then € <^ 1 and we 
have only a phase shift. That is, all collisions are 
adiabatic if V 2X < AE21. If adjacent levels are mixed, 
we will have phase shifts produced by slow (adiabatic) 
collisions and inelastic transitions produced by fast 
(nonadiabatic) collisions. 

For a series of collisions which do not overlap in 
time (or for which this overlap is negligible) the classi- 
cal form is preserved and the net result is a phase shift 
7] and a damping e which are obtained by adding up 
the phase shifts and damping constants from each of 
the collisions. When two collisions overlap, quantum 
mechanical mixing effects occur and the classical 
form is destroyed. 

d. Overlapping Lines 

If the lines a>io and (o 2 o are not isolated we must 
consider all terms in eq (51): 

C(V(t), t) *[p x d\J[,+p 2 d 2() d^T' vl \ exp (-ioiot) 

+ [p2d2()T 22 -\-pid U) do 2 T 2l ] exp (— ia)2ot). (63) 

The p\T'n term gives the probability that a radiator 
initially in the state E\ will remain in that state after a 
time t; the p 2 T\ 2 term gives the probability that a radia- 
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tor initially in the state E 2 will make a radiationless 
transition to E\ during the time t; together these two 
terms describe the probability, at the time t, that a 
radiator will emit a photon with a frequency near coio. 
The remaining terms in eq (63) describe the probability 
of emission near (020 (we must consider these terms 
when the lines o>io and 0*20 overlap). 

When t < £0, the radiator is unperturbed hence 
T' U = T! Z2 =1 and F; 2 = 7V=0; in this case, C(V(t), t) 
<* [pitffo exp (— icoiot) + pzdfco exp (— ia^ot)]. When the 
collision occurs, we get a probability for switching 
between the modes coio and cl^o as discussed earlier; 
that is, T[ 2 and T' n increase in amplitude while T' n and 
7^22 decrease. To show this effect, we will calculate the 
T' matrix elements for t > {to + r). 

The condition for overlapping lines is {AE 2 \/h) 
~Awi/2; since we may identify Aa)i/ 2 with v c where 
(Ijvc) is the mean free time between collisions and 
since r is usually much less than \jv c (at least for 
strong collisions), we must have (^E 2 \TJh)< 1. Noting 
the results of eqs (57) and (58), we see that when 
(A2?2iT/ft) < 1, the S matrix elements may easily be 
evaluated to all orders with the following results: 



Sn=S22 = cos (V^r/h) 
S\ 2 = S 2i =— i sin (V-nTJli). 



(64) 



Using eqs (55), (63), and (64) we may give an explicit 
expression for C{V(t), t) in terms of (V^rlh). To 
simplify the following discussion we will use C(V(t), t) 
= C x {V{t), t)+C 2 (V(t), 0, where 

Ci(V(t), t) ^p^tffoexp (— iojiot) cos (V v >Tlh) 

— idioidwl d 2 o) exp (— ia) 2 ot -\- ico 2 it ) sin (V 2i rlh)], 

(65) 

and C 2 (V(t), t) contains the two p 2 terms in eq (63). 
(C 2 may be obtained by interchanging the indices 2 
and 1 in C\.) We will discuss only C\ since the inter- 
pretation of C> is obviously identical. 

Ci(V(t), t) describes the behavior of a radiator 
which was initially in the state E\. Before the collision 



Ci(V(t),t) ^piJfoexp (-io) U) t) 



(66) 



and after the collision, t> £o + r, C\ is given by eq (65). 
cos 2 (V 2 irlh) gives the probability that the radiator 
will be in the state E\ after the collision, and sin 2 
{V 2 \rjfi) gives the probability that the collision will 
induce a radiationless transition to the state E 2 . If 
(V v >rlfi) = (77/2), the collision will definitely cause a 
transition to E 2 and we would have 

C x {V{t), t) oc pi^lo(rfio/^2o) exp (— ia) 2 ot+ io) 2 ito) 

(67) 

after the collision. Equation (67) describes a light 
train for the line 0*20 with a phase shift a> 2 i £0 just 
as in the classical oscillator discussed earlier (the 
factor diold 2 o represents the arbitrary initial phase 



difference between the two modes as well as an ampli- 
tude difference if one should exist). This result is 
identical with the result for a two mode classical oscil- 
lator because, when (Vi 2 r/h) equals some odd multiple 
of (7r/2), there is definitely a change of mode; for 
other values of (Vi 2 r/h\ the quantum mechanical 
expression, eq (65), shows only a probability for this 
mode switching effect. Notice that the mode switching 
or phase-memory effects are described by the off- 
diagonal 7" matrix elements. These off-diagonal 
elements were not needed in the description of iso- 
lated lines because the average phase shift o>2i£o 
is greater than unity (recall that the average value 
of Jo is the mean free time (l/^ c ) which is estimated 
by l/Acoi/2) and these elements vanish after the to 
average. 

In the case of overlapping lines, C{V(t), t) does 
not have the simple classical form which was obtained 
for isolated lines. In this case we would average 
C(V(t) y t) over collision parameters by defining a 
quantity e={V\ 2 TJh) and a function P(e) which gives 
the probability that the collision will produce a 
(V12T I ft) whose magnitude is e. The average over 
C(V(i), t) contains terms like 

I P(e) cos ede and I P(e) sin ede. 

If P{e) does not vary rapidly then we may say that 
collisions for which e > 1 effectively terminate the 
wave train (this argument is analogous to the treat- 
ment of large phase shifts in the classical theories). 



7.3. Strong Collisions 

In discussing the approximations relevant to the 
density matrix and the wave functions, our general 
procedure will be to show that each approximation 
breaks down only during a "strong collision." It 
is thus convenient at this point to clarify our meaning 
of strong collisions and to discuss briefly the dynamics 
of certain types of collisions. 

In the previous section it was shown that each col- 
lision produces a phase shift and changes the ampli- 
tude of the light train. In this manner, a series of 
collisions will eventually destroy the correlation 
between d(t) and d(0) and reduce C(t) to zero. It 
may happen that a single collision is capable of 
reducing C(t) to zero; such an event is defined to be 
a strong collision. From the results of the previous 
section we see that a strong collision is one which 
produces a large phase shift, rj > 1, or a large decrease 
in amplitude, e > 1. Such collisions terminate the 
light train and therefore satisfy the hypotheses of 
the Lorentz interruption theory when e > 1, or Weiss- 
kopfs phase shift theory when 77 > 1. The general 
program of classical path theories is to describe all 
but the strong collisions assuming that these will 
be treated by some other method such as a Lorentz- 
Weisskopf approximation [21]. 

From the results of the previous section, we see that 
a strong collision will occur when 
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(Vrlh)>\ 



(68) 



unless the level spacing is greater than V (cf eq. (61) 
using A£ 2 i > V n )- If AA21 > P21 we must replace the 
V in eq (68) by V, n = FjJAE 2l ; in this case, V ett is the 
actual perturbation of the radiator since we have used 
V u = 0. It is in this sense that we use eq (68) as a work- 
ing definition of strong collisions. It should be noted 
that our definition of a strong collision as one which 
terminates the light train is basically equivalent to 
the strong collisions in the impact theory [21 J which 
are said to occur when the terms in the S-matrix 
expansion (see eq (54)) are on the order of unity or 
larger. 

Most classical potential functions are inversely 
proportional to some power of the perturber position 



V(t)e*fh + 



JO 



)dt', where p,> is the radiator-per- 



turber separation at the time of closest approach 
and \{t) is the instantaneous perturber velocity. 
The maximum value of the potential where \r(t)\ 
is smallest, namely when r(0) = po. The half maximum 



occurs w 



hen 



FT 



= |Pd + 



Jo 



)dt'\=^ 2p () (which 



may be taken as a definition of the collision time r). 
lf\(t) changes very little during a collision we may use 
v(r) a* v(0) = v o or \r{r)\ =* |flo + TVo| =* 2p . Since 
Po* v = 0, the collision time is on the order of 

T=polvo (69) 

where vo is the velocity at the time of closest approach. 
It can be shown by a rather lengthy analysis of hyper- 
bolic trajectories, that this estimate is valid even when 
the radiator is an ion which accelerates a charged 
perturber. 

If the radiator is electrically neutral, we may esti- 
mate vo by the thermal velocity v av = V ikT/m. If both 
radiator and perturber are charged particles, we may 
estimate v () by 



vl=vl y — 2Z r Z t) e 2 lmp {) 



(70) 



where Z, and Z p denote the charge of the radiator and 
perturber respectively. 

8. Validity Criteria 

8.1. General 

An expression for the "thermal average" has been 
derived, eq (47), using a classical path treatment. 
This expression, or one like it, is the starting point 
for most of the modern semiclassical line broadening 
theories and it may be regarded as the end result of 
the classical path approximation. It is therefore worth- 
while to examine the validity criteria for this expression. 

It has been shown that the classical path approxi- 
mation, as it is used in line broadening theories, 
actually embodies several types of approximations. 
As we have stated them, they may be roughly divided 



into three classes; (1) the representation of the gas by 
a collection of noninteracting cells, (2) simplifications 
of the density matrix, and (3) a classical treatment of 
the perturber wave functions. 

In order to justify the mathematical model which 
is used by most line broadening theories, we have 
assumed that the average spacing between radiators 
is large enough that all direct and indirect interactions 
between radiators are negligible. (By indirect inter- 
actions we mean that one radiator may influence a 
perturber which is interacting with another radiator.) 
It is probably possible to justify the mathematical 
model with more relaxed assumptions, however, 
since our assumptions are satisfied for virtually all 
problems of interest, this may not be necessary. For 
resonance broadening in neutral gases, this approxi- 
mation may have to be reconsidered. 

It should be noted that the radiator Hamiltonian 
H a should describe translational motion of the radiators 
as well as their internal structure (i.e., bound states). 
In practice, however, one usually neglects this trans- 
lational motion and calculates only the discrete spec- 
trum. This procedure neglects Doppler shifts which 
also broaden a given line. The Doppler broadening 
may be reinstated, in an approximate manner, by 
means of convolution integrals [6J. This folding pro- 
cedure should be valid if the radiators trajectory, 
during the time of interest, is not appreciably altered 
by collisions with other particles. For most plasma 
line broadening problems this criterion is well satisfied 
however, for some foreign gas broadening problems 
one may have to treat Doppler effects more carefully. 

8.2. Density Matrix 

The approximations relevant to the density matrix 
are for the purpose of achieving a product form, eq 
(8), which permits us to factor the trace operations, 
as in eq (11), and thereby perform the average over 
perturber states without affecting the sum over states 
of the radiator. In this weak coupling approximation, 
the subsystems are statistically independent and this 
is what is meant by the neglect of back reaction as it 
applies to the statistics (density matrix). This approxi- 
mation has been made by replacing the term exp 
{—VjkT} in the Boltzmann factor by unity; that is, 
we have assumed that all V matrix elements (or at 
least those which pertain to the spectral line of in- 
terest) will satisfy V < AT. It should be noted that in 
obtaining this inequality we have tacitly assumed that 
the radiators and perturbers are in thermal equilibrium 
at a temperature T. For many line broadening experi- 
ments, however, the radiators and perturbers have 
different temperatures. Nonetheless, if the perturbers 
are in thermal equilibrium at a temperature 7)„ the 
weak coupling approximation will still be valid for 
V <^ kT p and the perturber statistics will be described 
by a Boltzmann factor as in eq (47). In such a case, 
the statistical weights for the atomic states (matrix 
elements of p {a) ) would be determined from the rate 
equations which describe the population and depopu- 
lation of atomic levels. We will henceforth take 
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V<kT, where T is the perturber temperature, as 
the validity criterion for the weak coupling approxi- 
mation. It is obvious that this inequality will be vio- 
lated for some interactions, however, it will be shown 
that the errors made by this approximation are impor- 
tant only in the far wings of a line where hAa) > kT. 

If the perturbation V> kT is produced by a static 
ion (one which does not move appreciably during the 
time of interest, see eq (52)) in a plasma, the Stark 
splitting is on the order of kT, and it is immediately 
obvious that this will influence only the far wings, 
Aco > kTjh, of a line. 

When the perturbation V > kT is produced during 
a collision with a dynamic perturber (we will con- 
sider electrons as representatives of dynamic per- 
turbers), the radiator-pert urber separation must be 
on the order of -K- dv =h/ VSmkT (the thermal wave- 
length, or average de Broglie wavelength, for the 
perturber) or smaller [22]. Using V=kT, po=^av 
and vo = v ax in eq (68), we see that, for neutral radia- 
tors, V > kT can occur only during a strong collision 
(for most neutral radiators the energy spacing between 
adjacent states is less than A^and it is not necessary 
to consider V eff as discussed in sec. 7.3). If a charged 
perturber collides with an oppositely charged radiating 
ion, it will be accelerated by the Coulomb attraction 



Fo, hence fo>^av. Using V=kT, po =y X a 



and 



Vo = v ay (volv aLV ) in eq (68), we find that (Vr/h) is on 
the order of (v a Jvo\ which is less than 1. Thus, 
in this case, V>kT need not correspond to a strong 
collision. In a collision involving a radiating ion and 
a charged perturber of like sign, the perturber will 
be decelerated, and it is unlikely that it will come 
close enough to the radiator to produce V>kT. 
To show this, we note that the smallest value of pc 
gi ven by eq (70) is (2e 2 /mpo)=v 2 av or po = X, av (2e 2 //i) 
vm/SkT; in order that po be on the order of ^X av , it 
is necessary that T be larger than 10 8 deg K. 

Regarding C{t) as a weighted sum of dd(£) terms 
as discussed in section 7.1, it is clear that the approxi- 
mate density matrix p i(l) p (l)) will incorrectly weight 
only those d-d(j) terms which correspond to a 
radiator initially undergoing a V<LkT collision (i.e., 
the density matrix weights each term according to the 
state at £ = 0). This does not preclude the possibility 
that V> kT collisions may occur at later times, t > 0, 
in the remaining, correctly weighted, d-d{t) terms. 
Such events will indeed occur, however, they are 
determined by the dynamics of the system and this 
has no connection with the statistical weighting (den- 
sity matrix) of the corresponding initial states. If we 
define a collision frequency vut for the V >kT collisions 
then the probability that such a collision will occur 
within a time r of £ = is Vkfr. The number of incor- 
rectly weighted d*d(0 terms will therefore be vuttNt, 
where Nr is the total number of terms. The number of 
correctly weighted terms will be (1 — Victt)Nt, and the 
error resulting from the use of the weak coupling den- 
sity matrix will be on the order of Vktrl (1 — Vkfr) * 
Since the duration, t, of a V > kT collision will be 
much smaller than the time between such collisions, 
we may replace (1 — Vkfr) by 1 and the error estimate 



becomes v^t. The error will be largest when all col- 
lisions produce V > kT; in this case we may estimate 
the collision frequency by Acoi/2, the halfwidth of the 
observed line. The largest value of r=po/i>o for a 
V^kT collision occurs for po = X av and Vo = v SLy (de- 
celerated perturbers, Vo < v av , are not likely to pro- 
duce V^ kT). The largest value of r is hjkT hence 
the largest error is on the order of hAcoi/z/kT. Since 
^Acoi/2 is much less than kT for virtually all lines of 
interest, this error is negligible. A more detailed 
analysis of the incorrectly weighted terms shows that 
this error estimate breaks down when Aco > Acoi/ 2 . In 
this case the time of interest is less than the average 
time between collisions (recall that Acoi/2 = VkT for 
the worst case), and some of the correctly weighted 
d-d(t) terms will not suffer a collision during the 
time of interest. These "unperturbed" terms will have 
little influence on the line shape in the region Aco> Acoi/ 2 
(their dominant influence is in the line center when 
Aco < Acoi/2). We must therefore compare the number 
of incorrectly weighted terms with (vkrl Aco) (1 — victt)Nt 
(the number of correctly weighted terms which have 
undergone a collision during the time of interest, 1/Aco). 
In this case the error estimate becomes hAco/kT which 
is small if Aco <^ kTjh. It must be emphasized that these 
error estimates are only an upper bound for the true 
error which, in most cases, is much less than 
(«Ao>W*T) or (hAio/kT). 

For neutral radiators V^kT occurs only during 
a strong collision. Since a strong collision completely 
disrupts the light train, the weighted sum of all 
d -d(t) terms which correspond to a strong collision 
at t = will have fallen to zero by the time these colli- 
sions are completed. That is, the halfwidth for the 
sum of all incorrectly weighted d • d(t) terms will 
be on the order of T = h/kT or smaller. The corre- 
sponding Fourier transform will thus be a constant for 
hAo) < kT. Since the incorrectly weighted terms 
add only a constant to 1(a)) for hAa) < kT, the approxi- 
mate density matrix will not affect the shape of 
/(co) for neutral radiators. Since most V> kT collisions 
are strong collisions even when the radiator is an ion, 
the above argument may be loosely applied to ion 
lines as well. We may therefore give the condition 
for validity of the weak coupling density matrix as 



hAa> <€ kT. 



(71) 



This condition was also obtained by Baranger [5,23] 
within the framework of a quantum mechanical impact 
theory. We note, in passing, that the region of validity 
for the impact approximation is given ToAco <^ 1 
where To is some representative collision time (see 
sec. 5.1 of ref. [4]). Since To > hjkT, we may state 
that when the impact approximation is valid the weak 
coupling density matrix may be used; the converse 
of this statement is not necessarily true. 

8.3. Wave Functions 

The dynamics of the system are needed only to 
calculate d(t) from d(0) for each given initial state; 
this requires a knowledge of T(t) but not ty (R, x, t). 
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That is, since we are working with a trace, eq (8), 
it is not necessary to know the wave functions 
M'(R, x, t) for the system (although it is necessary 
to define a set of basis functions). Semiclassical 
approximations have been made in the wave functions 
simply to permit the use of a classical potential 
function, F (/) (R, x (?) , v (/) , t), in the time development 
operator. These semiclassical approximations, which 
are the hallmark of the classical path approximation, 
permit a simple physical interpretation of the dynamics 
of the system. More powerful quantum mechanical 
treatments are available [3,4]; however, in many 
cases, the results of these jheories differ only slightly 
from the classical path treatment and the physical in- 
sight afforded by the latter is frequently advantageous. 
The wave packet assumption was used to obtain 
the relation 



K(R, xV (i) (x, t)= F (/) (R, x<*>, v (/) , f)<p (/) (x, r), (72) 

where V {i) is a classical function of the classical 
position and velocity coordinates x (0 and v 0) . V(R, x) 
can be expressed as a sum of interactions F(R, Xj) 
between the radiator and each perturber (recall that 
x=(xi, X2, . . . x.v)), thus we need only require that 



F(R, XjV'Hxj, 0=P (i) (R, x«> v<| t)<p«Xxj,t) (73) 

when F(R, xj) 7^ in order to obtain eq (72). Mes- 
siah [24] has shown that eq (73) is obtained if the 
perturber is represented by a classical wave packet 
[25, 26] and if the potential field seen by the per- 
turber does not vary appreciably (as compared with 
the perturbers kinetic energy) over the width of the 
wave packet. 

The conditions for validity of eq (73) are discussed 
in detail by Messiah and we shall only outline this 
discussion as it applies to our problem. We must 
show; (1) that it is possible to construct a wave packet 
which does not overlap with the wave functions for 
the other particles, (2) that this wave packet does not 
spread appreciably during a collision (i.e., when 
F(R, Xj) is nonzero), and (3) that the potential seen 
by the perturber does not vary appreciably over the 
width of a wave packet. The potential seen by the 
perturbers may be obtained by iterating the Schro- 
dinger equations for x (/) an( l <£ (/) to the next higher 
order as discussed in section 5. This potential will 
be the matrix elements of F(R, x) between atomic 
states or, in some cases, it will be V e{{ as discussed 
in section 7.3, (or some other type of optical potential) 
[27]. The momentum of a perturber moving in the 

potential F(R, x) is given by p=V2m(E—V), 
where E denotes the perturber energy. From this 
v expression we see that Ap = mAV/p, where Ap 
and AV denote the variation in p and Fover the width 
of a wave packet. The condition that AV < p 2 /2m 
is thus equivalent to 



which is the usual condition for validity of WKB calcu- 
lations (see eq (7.13) of ref. [28] noting that k(dp/dx) 
may be replaced by Ap). We will assume for the mo- 
ment that this inequality can be satisfied, and it will 
be shown that the other validity conditions are easily 
met. 

We wish to represent those perturbers which are 
not interacting with the radiator by Gaussian wave- 
packets of width Ax f , and momentum spread Ap^ <^ p av 
at some instant of time. For this wave packet construc- 
tion we must require that Ax g be much less than the 
average particle spacing n~ 113 (n is the perturber 
density) so that the wave packets do not overlap. 
Using AxgApg — h (for a Gaussian wavepacket), we see 
that the condition Ap,, < p av also requires Ax g > h/p ay 
=A av . We must therefore be able to satisfy. 



-ft. AV <^ Ax (l < // 



(75) 



Ap <$ p 



(74) 



For virtually all problems of interest we have -ft av < 
n~ 113 , hence it will be possible to construct Gaussian 
wave packets for the overwhelming majority of per- 
turbers at some instant of time. 

Since static perturbers (e.g., ions in a plasma) are 
regarded as stationary during the time of interest, 
it is only necessary that they be represented by wave 
packets at some instant of time. As discussed above, 
this will be possible for the overwhelming majority 
of perturbers if ^ flr <^ n~ 113 . The treatment of static 
perturbers will be discussed further in II and for the 
remainder of this section we will consider only dynamic 
perturbers. 

During a collision a wave packet will generally in- 
crease both its width, Ax, and its momentum spread, 
Ap. If Ap < p (or AV < fr/2m n which will be verified 
shortly), the increase in momentum spread will be 
negligible, and for a collision of duration T=p /i;o, 
the wave packet will increase its width to a maximum 
of Ax T = (Ajc-frAp/m). To prevent overlap with the 
radiators wave function during the collision Ax T 
should satisfy (Ajt T + /? ) < p , where R denotes the 
effective "extent" of the perturber wave function. 
Using T=po/vo and Ap<^p, we see thai (rAp/m)<^po 
hence we need onlv require ( A # -+-/?<>) < p<>. 

For a collision involving a neutral radiator the 
perturber momentum may be taken to be p av , hence 
we may use Ax — Ax ( , and Ap — Ap () (which will then 
satisfy Ap <^ p). Since Ax may be within an order of 
magnitude of ft dv (and still satisfy Ax>K av ), we may 
obtain an order of magnitude estimate of the overlap 
condition, (Ax + R<j) < p , by considering (*fc av + /? ) < po- 
If R <^ av , the wave functions will overlap only when 
po<ft a v> an d f° r neutral radiators, this produces a 
strong collision (see sec. 8.2). If R () >1{ ax , a more 
detailed analysis of the perturbation potential is 
required; as an example we will consider the Stark 
broadening interaction, V—e 2 Rolp%, for which 
(Vr/h) - (Rolp»)(e 2 lhv av ). Since (e 2 lhv av ) > 1 for 
temperatures of 10 6 K or less (for electron perturbers), 
we see that the overlap condition, R {) > p<>, is violated 
only during a strong collision. Thus, for neutral radia- 
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tors, the wave functions will overlap only during a strong 
collision. In the case where a radiating ion accelerates a 
charged perturber the deBroglie wavelength of the 
perturber is decreased, &o = hlmvo<^ av . If vo is not 
much larger than i? av , there is essentially no difference 
between this case and the neutral radiator discussed 
above. If z;o^>z; av the situation is somewhat different. 
In this case the perturber will come very close to the 
radiator and it is necessary to construct the initial 
wave packet with a width Ax much smaller than ft av . 
Since -X.o^'^av^ we construct the wave packet so that 
Ajc >Xo which, using A;tAp > ft, implies that Ap < p () . 
Such a wave packet will not hold together very long, 
however, it will hold together for a time on the order 
of t= po/vo which is all that is required by eq (73). 
In this case we will have "Xo^ Ro and the wave func- 
tions will overlap when po < Ro. The perturbation 
potential will be 

V^ZeWo/pl and (Vr/h) = (Z<?Rolhp v<>); 

we do not need to consider a V e{{ as discussed in 
section 7.3 because when p () < /? , V is comparable to 
the ionization energy for ions [29]. If vo > v ay then, 
from eq (70), we have vo — (2Ze 2 /mp ) 1/2 and 

(Ft/*) - [(Z/2) (fl 2 /po«,) ] 1/2 

where ao= h 2 jme 2 is the Bohr radius (we have chosen m 
as the electron mass; for an ion perturber (Vr/h) 
would be even larger). Since Ro is larger than a , we 
see that the overlap of wave functions again produces 
a strong collision. 

It is interesting to note that, for Stark broadening 
in plasmas we have {Vr/h) — (Ze 2 RolhpoVo)> (ZRolaot), 
where /= mpoVo/h is the angular momentum of the 
perturber. Since Ro > ao, we will have strong colli- 
sions when / <; 1. Other interaction potentials may be 
treated in a similar manner with the same general 
results. In any case, the wave functions will not overlap 
if %) < p ; using ho = h/mvo and l = mpoVo/h, we find 
that the radiator and perturber will be distinct (i.e., 
their wave functions do not overlap) when 



/>!. 



(76) 



It is also interesting to note that the perturbers will 
look like classical particles as they pass the radiator 
only if Ax < x and Ap <^ p; using A^cAp — ft, we again 
obtain the inequality h < mpoVo or / > 1. We note 
further that A/ = ± 1 in a dipole transition, thus A/ < / 
and it is in this sense that the angular momentum of a 
perturber is conserved during a collision, as it must be 
if we are to use classical mechanics to describe the 
perturbers. 

We have now shown that the validity criteria for 
the use of classical wavepackets are all easily satisfied 
if Ap <^ p, or equivalently, AV < p 2 /2m. That is, the 
essential validity criterion for the classical wave- 
packets is the same as the condition for validity of a 
WKB calculation. This was also noted by Baranger [4] 
who found that the results of a quantum mechanical 



impact theory reduce to the classical path results 
when a WKB treatment of the wave functions is 
employed. 

To show that the condition AV < p-\2m is satisfied 
we note that V increases monotonically as the radiator- 
perturber separation is decreased, thus AV ^ V and it 
is sufficient to show that V < p 2 /2m. Again V may be 
replaced by F eff if V is not strong enough to mix 
adjacent states. Using p 2 /2m = mv-J2, T = po/t>o, and 
l=mpoVolfi, we obtain 



(Vr/hl) <1. 



(77) 



Noting eqs (68) and (76), it is clear that this inequality 
can be violated only during a strong collision (recall 
that / < 1 produces a strong collision). 

Equation (73) requires that a perturber be repre- 
sented by a wave packet when it collides with the 
radiator. If we let (p(*> (xj) denote a Gaussian wave 
packet for the yth particle then the wave function 
(p^fxj, t)=Tj{t)<p { '\xj y 0), where Tj is the time devel- 
opment operator for this particle, will be a Gaussian 
packet at some time tj if we construct c^ (/) (xj, 0) such 
that 

<p«>( Xj , 0)=r j (-^ ( i ) .p.(x j ) (note that T^j)T^-tj)=l). 

In this manner we can arrange the wave function 
<£> (,) (x, t) = n#> (,) (xj, t) in such a way that all particles 
which collide with the radiator during the time of 
interest are represented by Gaussian wavepackets 
at the time of their collision. In general these particles 
will not be represented by wavepackets before or 
after their collisions, but this does not matter as far 
as eq (73) is concerned. 

We have now shown that it is possible to represent 
the perturbers by. classical wavepackets as required 
by eq (73) unless a strong collision occurs (in which 
case we do not use wavepackets). It is interesting to 
note that it may be possible to treat some strong col- 
lisions by classical path methods. That is, eq (77) could 
be satisfied when (Vr/h) ^ 1 if I is sufficiently large. 
If one wished to do this, it would be necessary to 
reexamine the validity criteria in more detail. For the 
purposes of this paper, however, it is sufficient to 
state that the use of classical wavepackets is justified 
for all collisions which are not strong collisions. 



8.4. An Approximate Classical Potential 

Although we have justified the use of a classical 
potential as given by eq (73), we have not yet specified 
its functional form. If back reaction effects in <p® 
are negligible, the potential V^ will result from 
a perturber passing by the radiator on a straight 
line trajectory (or a hyperbola if both radiator and 
perturber are charged particles). If back reaction 
effects are important, the trajectory will be a rather 
complicated function of velocity, acceleration, etc. 
A wavepacket formed by WKB calculations (the 
iterative procedure discussed in sec. 5) may be 
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able to account for the true classical trajectory, how- 
ever, the perturber wave function (p (i > would no 
longer be an //,, eigenfunction. While this may not 
matter for some scattering problems, it makes a 
great deal of difference in line broadening. 

In section 6 we used the fact that <p^ is an //,, 
eigenfunction to obtain the classical energy expression 
£^(x (i) ,v (<) ) in eq (44). If the perturber wave function 
were not an energy eigenfunction, the perturber 
density matrix would not be "diagonal in wave packet 
states" and the average over perturber states, eq (41), 
would be intractable in a wave packet representation. 
We must therefore show: (1) that <p (/) may be an 
Hj, eigenfunction (i.e., back reaction effects in <p® 
are negligible), and (2) that the energy eigenvalues 
may be represented by £W(x*'\ v (l) ) as stated in 
eq (44). 

Since energy eigenvalues are not time dependent, 
we need only show that E$ represents the H v eigen- 
value at some instant of time. This problem is essen- 
tially the same as eq (73) except that all perturbers 
must be represented by classical wave packets at 
some instant. Since we have a large number of per- 
turbers in our system (a cell), there will always be 
some wave functions which overlap; however, for 
virtually all problems of interest we have ^C av < //~ 1/: * 
and it is possible to construct wave packets for the 
overwhelming majority of perturbers at some instant 
of time (cf. eq (75)). It is in this sense that eq (44) 
is satisfied. 

To show that back reaction effects in <p^ are negli- 
gible, we consider the perturbation F* (x 0) , \ (l) ,t) 
which results when a perturber moves past a neutral 
radiator on a straight line trajectory (the trajectory 
will be a hyperbola if both perturber and radiator are 
charged particles). The true perturber trajectory 
will deviate from a straight line (or hyperbola) because 
of accelerations (back reaction) resulting from the 
potential V seen by the perturbers (this potential 
may be V^ as discussed in sec. 7.3). The approximate 
trajectory will be valid if this deviation does not appre- 
ciably alter the potential V® seen by the radiator 
during the collision (i.e., during the time r — po/fo, 
or the time of interest 1/Aco, if this should be shorter). 
V {i) is influenced by the perturber trajectory through 
its dependence on r (t), the position of the perturber 
at the time t. If the deviation, Ar(£), produced by 
the perturbation is small enough that r(t) ■+■ Ar(t) — r(t), 
the approximate trajectory will provide a good estimate 
of the true perturbation potential (whether the radiator 
is an ion or a neutral). 

For a time on the order of r, the magnitude of 
the deviation may be estimated by Ar = |F | r 2 /2m 
= \Vy\7 2 /m where ¥= — W is the back reaction 
force which causes the perturber to deviate from the 
approximate trajectory. Since V is inversely propor- 
tional to some power of r(t) we have \VV\ — V/r{t). 
Using r(t) ^ p<> we have 



where r=po/Fo, l=mp vo/h, and / > 1. This result 
shows that (Ar/r)<^ 1 unless a strong collision occurs. 
The above argument shows that the approximate 
trajectory is valid (unless a strong collision occurs) 
for elastic collisions. If an inelastic collision should 
occur, the radiator and perturber will exchange an 
amount of energy AE. This will produce a sudden 
change in the perturber velocity, AE = mvAv^ which 
will give rise to a deviation Ar ^ tAv in the trajectory. 
Again using To = po/t>o, l = mvopolh, and /> 1, we obtain 



Ar Av _ AE rAE tAE 
r ~~ v { ) imvf) hi h 



(79) 



Ar Vt 2 Vt Vt 

^ < <^ 

/• 2 m pi hi h 



(78) 



In section 7.2 it was shown that the probability for 
inelastic transitions is essentially zero unless 
{rAE/h) < 1 (for a more detailed treatment see ref. 
[30]). .Thus the approximate trajectory is valid for both 
elastic and inelastic collisions. 

It is interesting to note, that the change in the per- 
turbers kinetic energy during the collision is a small 
fraction of its total kinetic energy. 

9. Summary and Comments 

It has been found that if the approximations which 
are used in the classical path treatment of line broaden- 
ing break down, then a strong collision, defined by eq 
(68), is occuring. It should be noted that the converse 
of this statement is not necessarily true; that is, it 
may be possible to extend the classical path treat- 
ment of some problems into the region where [Vrlh) — 1. 
In this case the validity conditions must be considered 
in more detail. 

In the application of classical path methods, strong 
collisions are usually treated by interruption approxi- 
mations based on the Lorentz-Weisskopf theory; the 
influence of all other collisions is evaluated by means 
of the time development operator T\l\t) as discussed 
in section 6. 

In section 8.4 it was shown that back reaction effects 
on the perturber trajectory are negligible whenever 
a classical path calculation is justified. That is, the 
validity criteria for the use of classical wave packets 
break down at the same time as (or in some cases 
before) the validity conditions for the use of unper- 
turbed trajectories. 

It was noted that our definition of strong collisions, 
eq (68), is essentially the same as the result obtained 
by setting the lead terms in an S-matrix expansion to 
unity. However, using this S-matrix expansion until 
the lead terms are on the order of unity may be a very 
poor approximation of the S-matrix; indeed, most errors 
attributed to strong collisions stem from this series 
approximation. To avoid possible confusion, we wish 
to emphasize that this approximation concerning the 
S-matrix is not part of the classical path approxima- 
tion and any errors resulting from it have no bearing 
on the validity of classical path methods. It is therefore 
quite possible to improve upon these so-called strong 
collision problems within the framework of the classi- 
cal path approach. The series approximation has been 
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improved for Lyman alpha lines [31], by summing 
terms in the series and, for hydrogen lines in general 
[32], by an approximate treatment of the exponential 
definition of the S-matrix. 

The region of validity for the weak coupling density 
matrix has been given as Aco < kTjh. We wish to em- 
phasize once again that the approximations relevant 
to this density matrix are not influenced by the semi- 
classical treatment of the wave functions which is a 
part of the classical path approximation. The weak 
coupling density matrix may be (and has been [3, 4]) 
used in more rigorous quantum mechanical theories. 
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